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Abstract 

We  consider  the  problem  of  tracking  multiple  maneuvering  targets  in  the  presence  of  clutter 
using  switching  multiple  target  motion  models.  A  novel  suboptimal  filtering  algorithm  is  devel¬ 
oped  by  applying  the  basic  interacting  multiple  model  (IMM)  approach  and  joint  probability 
data  association  technique.  But  unlike  the  standard  single  scan  joint  probabilistic  data  associa¬ 
tion  (JPDA)  approach,  we  exploit  a  multiscan  joint  probabilistic  data  association  (Mscan- JPDA) 
approach  to  solve  the  data  association  problem.  The  algorithm  is  illustrated  via  a  simulation 
example  involving  tracking  of  three  maneuvering  targets  and  a  multiscan  data  window  of  length 
two. 

Keywords:  Multitarget  tracking;  interacting  multiple  model  (IMM)  algorithm;  multiscan  joint 
probabilistic  data  association;  state  estimation. 

I  Introduction 

We  consider  the  problem  of  tracking  multiple  maneuvering  targets  in  presence  of  clutter  using 
switching  multiple  target  motion  models.  This  class  of  problem  has  received  considerable  attention 
in  the  literature  [2], [5], [6], [8], [10], [23].  The  switching  multiple  model  approach  has  been  found  to 
be  very  effective  in  modeling  maneuvering  targets  [2]- [7], [9], [12].  In  this  approach  various  modes  of 
target  motion  are  represented  by  distinct  kinematic  models,  and  in  a  Bayesian  framework,  the  target 
maneuvers  are  modeled  by  switching  among  these  models  controlled  by  a  Markov  chain.  While 
tracking  multiple  targets  in  the  presence  of  clutter,  one  has  to  solve  the  problem  of  measurement 
origin  uncertainty,  i.e.  how  to  associate  the  data  available  at  the  sensor(s)  with  various  targets  or 
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clutter  (false  measurements).  This  problem  of  data  association  has  already  been  considered  and 
effectively  solved  in  the  Bayesian  framework  by  using  probabilistic  data  association  (PDA)  [1),[2],[5] 
in  the  case  of  a  single  target  and  joint  probabilistic  data  association  (JPDA)  [2], [5], [10]  in  the  case 
of  multiple  targets.  The  standard  JPDA  algorithm  uses  only  a  single  (latest)  scan  data  available  at 
the  sensors  and  the  state  of  the  target  is  updated  by  using  a  weighted  sum  of  measurements  which 
could  have  reasonably  originated  from  the  target  under  consideration.  To  use  more  information 
to  solve  the  data  association  problem,  the  idea  of  using  multiple  scans  of  data  (current  and  past 
scans)  seems  to  have  been  initially  proposed  by  Drummond  [19]- [21].  Drummond  describes  some 
practical  issues  involved  but  does  not  discuss  detailed  problem  formulation  and  technical  issues  for 
multiscan  JPDA.  Roecker  [14]  has  extended  Drummond’s  ideas  where  he  has  discussed  problem 
formulation  and  solution  in  some  detail.  In  a  simulation  example  presented  in  [14]  it  has  been 
shown  that  performance  improvement  in  multiple  target  tracking  can  be  achieved  via  multiscan 
JPDA  as  compared  to  single  scan  JPDA  with  most  of  the  improvement  gains  achieved  via  a  window 
size  of  2  or  3  scans. 

In  [14]  only  non-maneuvering  targets  (i.e.  one  model  per  target)  have  been  considered.  In  this 
paper,  we  extend  Roecker’s  approach  to  highly  maneuvering  targets  where  we  allow  multiple  kine¬ 
matic  motion  models  per  target.  A  novel  suboptimal  filtering  algorithm  is  developed  by  applying 
the  basic  interacting  multiple  model  (IMM)  approach  and  multiple  scan  joint  probability  data  asso¬ 
ciation  (Mscan-JPDA)  technique.  The  algorithm  is  illustrated  via  a  simulation  example  involving 
tracking  of  three  maneuvering  targets  and  multiscan  data  window  of  length  two. 

The  paper  is  organized  as  follows.  The  basic  multiscan  JPDA  problem  is  explained  in  Sec. 
II  followed  by  the  presentation  of  the  problem  formulation  in  Sec.  III.  The  proposed  IMM-based 
multiscan  JPDA  algorithm  is  described  in  Sec.  IV  for  the  case  of  a  sliding  raultiscan  window  of  size 
2.  A  computer  simulation  example  is  presented  in  Sec.  V. 

II  Multiscan  Joint  Probabilistic  Data  Association 

The  first  data  association  problem  for  target  tracking  was  addressed  by  Bar-Shalom  and  Tse  [1], 
in  which  they  considered  single  target  tracking  in  the  presence  of  clutter.  They  developed  the 
PDA,  a  method  of  data  association,  in  which  a  probability  weight  for  each  measurement-to-target 
association  is  computed  and  the  state  of  the  target  is  updated  with  a  weighted  sum  of  measure¬ 
ments.  For  multiple  target  tracking  in  presence  of  clutter,  the  JPDA  algorithm  has  been  developed 
[2], [5], [10].  Unlike  PDA,  JPDA  computes  the  probabilistic  weight  for  measurement-to-target  as¬ 
sociation  jointly  across  the  set  of  all  targets  and  clutter.  Basically  it  defines  all  the  feasible  joint 
events  for  the  known  number  of  targets  and  clutter.  Each  feasible  joint  event  is  a  unique  event  that 
represents  the  association  of  measurements  to  targets  and  clutter.  The  probability  of  each  joint 
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event  is  evaluated  and  the  state  of  the  target  is  updated  with  a  weighted  sum  of  joint  events  in 
which  a  measurement  is  associated  with  the  target  under  consideration. 

A  disadvantage  of  JPDA  is  that  it  uses  only  the  data  present  in  the  current  scan.  The  idea 
of  multiscan  JPDA  (n-scan-back)  stems  from  the  fact  that  if  more  scans  of  data  are  used,  more 
information  would  be  available  for  data  association  and  hence  for  computing  the  probabilistic 
weights.  One  can  use  all  combinations  of  measurements  beginning  from  the  initial  time  to  the 
present  time,  so  that  all  information  can  be  utilized  to  compute  the  probabilistic  weights  and  hence 
to  produce  the  best  state  update.  This  “optimal  bayesian  filter”  has  been  considered  previously 
[2], [5].  But  practically  it  is  not  feasible  to  implement  due  to  memory  usage  and  exponentially 
increasing  computations.  Hence  practically  feasible  n-scan-back  or  multiscan  JPDA  approach  will 
be  considered  with  a  fixed  size  sliding  window  of  n  scans  (current  and  past  n  —  1  scans). 

The  basic  idea  of  multiScan  JPDA  can  be  explained  as  follows.  In  a  single  scan  JPDA  we  define 
single  scan  joint  events.  Similarly  in  the  multiscan  scenario  we  define  multiple-scan  joint  events 
as  follows.  A  marginal  association  event  6ir(k)  is  said  to  be  effective  at  time  scan  k  when  the 
validated  measurement  is  associated  with  (i.e.  originates  from)  target  r  (r  =  0, 1,  •  •  • ,  N  where 
r  =  0  means  that  the  measurement  is  caused  by  clutter).  Assuming  that  there  are  no  unresolved 
measurements,  a  joint  association  event  0*  is  said  to  be  effective  when  a  set  of  marginal  events 
(0,r(/c)}  holds  true  simultaneously.  That  is,  ©fc  =  f)£Li  6in(k)  where  r*  is  the  index  of  the  target 
to  which  measurement  y]^  is  associated  in  the  event  under  consideration,  (i  =  1, 2,  •  •  • ,  m).  In  the 
multiscan  case  with  a  scan  window  size  L  (A-scan-back)  and  ks  =  k  —  L  4-  s,  we  define  multiscan 
joint  events 

Lm 

®kL  =  Pi  P|  0iria(&s) 

5=1  i=l 

where  8iris(ks)  is  the  marginal  association  event  that  at  time  scan  ks,  ith  the  validated  measurement 
yj^  is  associated  with  target  r;s.  Let  |{0fc}|  denote  the  total  number  of  feasible  joint  events  in  the 
single  scan  case.  In  the  multiscan  case  we  get  total  number  of  multiscan  feasible  joint  events  as 
|{0fc}|  x  |{0fc_i}|  x  •  ••  x  |{©fc_£+i}|,  derived  from  a  Cartesian  product  of  joint  events  present  in 
all  scans  considered  in  the  scan  window.  As  one  can  see,  even  for  a  multiscan  window  of  length  two 
or  three,  the  number  of  feasible  multiscan  joint  events  grow  exponentially.  Rocker  [14]  and  Poore 
[18]  claim  that  a  multiscan  window  of  length  two  or  three  is  generally  enough;  any  further  increase 
in  the  window  length  achieves  only  marginal  performance  improvement. 

Once  we  define  the  feasible  multiscan  joint  events,  the  next  step  is  to  compute  the  probability 
weights  for  these  events.  The  heart  of  the  algorithm  is  finding  these  probability  weights  which  is 
discussed  in  detail  in  the  sequel. 
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Ill  Problem  Formulation 


Assume  that  there  are  total  N  targets  with  the  target  set  denoted  as  Tn  :=  {1, 2,  •  •  • ,  N}.  Assume 
that  the  dynamics  of  each  target  can  be  modeled  as  one  of  the  n  hypothesized  models.  The  model 
set  is  denoted  as  Mn  :=  {1,2,  •  •  •  ,n}.  For  target  r  ( r  G  Tn),  the  event  that  model  i  is  in  effect 
during  the  sampling  period  (t*_i,tfc]  will  be  denoted  by  Mxk{r).  Although  all  the  targets  share  a 
common  model  set,  any  two  targets  may  be  in  different  motion  status  from  time  to  time. 

For  the  j-th  hypothesized  model  (mode),  the  state  dynamics  and  measurements  of  target  r 
(r  G  Tn)  are  modeled  as 

Xk(r)  =  +  Gi_i(r)«fc-i(r)  (!) 

and 

Zk{r )  =  h3{xk{r))  +  w{(r)  (2) 

where  Xk(r)  is  the  system  state  of  target  r  at  tk  and  of  dimension  nx  (assuming  all  targets  share 
a  common  state  space),  zk(r)  is  the  (true)  measurement  vector  (i.e.  due  to  target  r)  at  tk  and  of 
dimension  nz,  F^_j(r)  and  G^._1(r)  are  the  system  matrices  when  model  j  is  in  effect  over  the 
sampling  period  {tk-i,tk]  for  target  r  and  /iJ  is  the  nonlinear  transformation  of  xk(r)  to  zk(r)  for 
model  j.  A  first-order  linearized  version  of  (2)  is  given  by 

zk{r)  =  H3k(r)xk(r)  +  toj(r)  (3) 


where  Hk(r)  is  the  Jacobian  matrix  of  evaluated  at  some  value  of  the  estimate  of  state  xk[r)  (see 
Sec.  IV).  The  nature  of  the  system  state,  the  various  matrices  in  (1)  and  (3),  and  the  measurements 
is  specified  in  more  detail  in  Sec.  V.  The  process  noise  v3k_1(r)  and  the  measurement  noise  vPk{r)  are 
mutually  uncorrelated  zero-mean  white  Gaussian  processes  with  covariance  matrices  Qk_1  (same  for 
all  targets)  and  R3k  (same  for  all  targets),  respectively.  At  the  initial  time  to,  the  initial  conditions 
for  the  system  state  of  target  r  under  each  model  j  are  assumed  to  be  Gaussian  random  variables 
with  the  known  mean  x30(r)  and  the  known  covariance  Pq  (r).  The  probability  of  target  r  in  model 
j  at  to,  Mo(r)  =  (r ) } ,  is  also  assumed  to  be  known.  The  switching  from  model  M^_j(r) 

to  model  Mk(r)  is  governed  by  a  finite-state  stationary  Markov  chain  (same  for  all  targets)  with 
known  transition  probabilities  pij  =  P{M3k(r)\Mtk_1(r)}.  Henceforth,  tk  will  be  denoted  by  k. 

The  following  notations  and  definitions  are  used  regarding  the  measurements.  Note  that,  in 
general,  at  any  time  k,  some  measurements  may  be  due  to  clutter  and  some  due  to  the  target,  i.e. 
there  can  be  more  than  a  single  measurement  at  time  k.  The  measurement  set  (not  yet  validated) 
generated  at  time  k  is  denoted  as 


Zk  :=  •  •  • ,  z ^ } 
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where  m  is  the  number  of  measurements  generated  at  time  k.  Variable  (i  =  1,  •  •  • ,  m)  is  the 
ith  measurement  within  the  set.  The  validated  set  of  measurements  at  time  k  will  be  denoted  by 
Yk,  containing  m  (<  m)  measurement  vectors.  The  cumulative  set  of  validated  measurements  up 
to  time  k  is  denoted  as 

Zk  =  {YuY2,---,Yk}.  (5) 

With  this  problem  formulation,  we  make  the  following  (standard)  assumptions  before  presenting 
the  detailed  algorithm. 

Assumptions: 

(1)  It  is  assumed  that  the  number  of  targets  ( N )  is  known  and  that  for  each  target  track  has 
been  initiated,  and  our  objective  is  to  maintain  the  tracks. 

(2)  Assuming  there  are  no  unresolved  measurements  (i.e.  measurements  associated  with  two  or 
more  targets  simultaneously),  any  measurement  therefore  is  either  associated  with  a  single 
target  or  caused  by  clutter. 

(3)  Clutter  is  modeled  as  independently  and  identically  distributed  (i.i.d.)  with  uniform  spatial 
distribution  over  the  entire  validation  region  (across  all  targets). 

(4)  State  estimate  of  individual  targets  conditioned  on  the  modes,  joint  events  and  set  of  mea¬ 
surements  are  mutually  independent  and  Gaussian  distributed  i.e.  states  of  the  targets  are 
not  coupled  and  estimation  is  carried  out  independently. 

(5)  Multiscan  window  of  length  two  will  be  used  to  compute  multiscan  joint  probabilities.  Ex¬ 
tension  to  higher  lengths  is  straightforward  but  tedious. 

The  goal  is  to  find  the  filtered  state  estimate  for  target  r  (r  €  Tn) 

xklk(r)=E{xk(r)\Zk}  (6) 

and  the  associated  error  covariance  matrix 

Pk\k(r )  =  E{[xk{r )  -  £fc|fc(r)][z*(r)  -  xk\k{r)]'\Zk}  (7) 

where  denotes  the  transpose  of  Xfc(r). 

IV  IMM/Mscan-JPDA  Filtering  Algorithm 

We  now  extend  the  single  scan  IMM/JPDA  filtering  algorithm  of  [8]  to  apply  to  the  multiscan 
case.  As  in  [14]  we  will  follow  a  sliding  window  multiscan  approach.  The  approach  of  [8],  in  turn, 
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is  based  on  the  approaches  of  [2],  [5],  [10],  [11].  As  the  IMM/PDAF  algorithm  is  well-explained 
in  [5,  Sec.  4.5]  and  [11],  the  JPDAF  algorithm  is  well-explained  in  [2,  Sec.  9.3]  and  [5,  Sec.  6.2], 
and  the  IMM/JPDA  filter  is  given  in  detail  in  [8]  (where  all  the  underlying  assumptions  and 
approximations  may  be  found  in  further  detail),  we  will  only  briefly  outline  the  basic  steps  in 
“one  cycle”  (i.e.  processing  needed  to  update  for  a  new  set  of  measurements  and  a  new  multiscan 
window)  of  the  IMM/JPDA  multiscan  filter.  We  assume  that  the  scan  window  size  is' two.  Given 
state  estimate  at  time  k  —  1  based  on  data  up  to  time  k  —  1,  in  Sec.  IV.l  we  provide  first  scan  steps 
(using  data  up  to  time  k)  and  in  Sec.  IV.2  we  provide  the  second  scan  steps  (using  data  up  to  time 
k  + 1). 

Assumed  available:  Given  the  state  estimate  $]._1|fc_1(r)  :=  E{xk-i{r)\Mk_l(r)>Zk~1},  the 
associated  covariance  i^_1|fc_1(r)  and  the  conditional  mode  probability  /4-i(r)  =  P[M/_1(r)|2fc_1] 
at  time  k  —  1  for  each  mode  j  G  Mn  and  each  target  r  6  Tn- 

IV.l  First  Scan  Steps 

Step  1.1.  Interaction  —  mixing  of  the  estimate  from  the  previous  time  (Vj  G  Mn, 

Vr  G  TN): 

n 

predicted  mode  probability:  :=  •F>{Affc(r)|.£fe-1}  =  ^,Pijtik-i(r)’  (8) 

1=1 

mixing  probability:  ^ 3{r )  :=  Zk~1}  =  VijiA-iir)/ n{~ (r).  (9) 

mixed  estimate:  :=  E{xk-X{r)\Ml{r),Zk~1}  =  (0/^(0-  (10) 

i= 1 

covariance  of  the  mixed  estimate: 

==  *{[**-iM  -  -  C,|t-iW]'|M’(r),^-'} 

=  +  K-iit-ifr)  -  ^Lin-i(r)][®*-ii»-i(r)  -  *°’-i|t-i(r)]  Vi,(r)-  (u) 

t=l 

Step  1.2.  Predicted  state  (Vj  G  Mn,  Vr  G  Tat): 

State  prediction:  4|fc-i (r)  E{xk(r)\M3k{r),  Zk~x }  =  ■  (12) 

State  prediction  error  covariance: 

Pk\k-i(r)  =  Ei[xk(r)  ~  ££|fc_i(r)]Mr)  -  ^_i(r)]'|Af/(r),2:fe-1} 

=  (i3) 
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Using  (2)  and  (12),  the  mode-conditioned  predicted  measurement  of  target  r  is 

*k(0  :=  (»*))•  (14) 

Using  the  linearized  version  (3),  the  covariance  of  the  mode-conditioned  residual  v^l\r)  :=  z^k  — 
z{(r)  is  given  by 

Sl(r)  :=  E{4U(r)viU(r)'\Mi(r),Zk-'}  =  (r)  +  R{  (15) 

where  H3k(r)  is  the  first  order  derivative  (Jacobian  matrix)  of  h3(.)  at  x(r).  Note  that  (15) 
assumes  that  originates  from  the  target  r. 

Step  1.3.  Measurement  validation:  There  are  two  steps  to  measurement  validation.  First 
perform  measurement  validation  for  each  target  r  (r  6  Tn)  separately.  For  target  r,  the  validation 
region  is  taken  to  be  the  same  for  all  models,  i.e.,  as  the  largest  of  them.  Let  (|A|  =  det(A)) 


Jr  :=  arg  j  max  \S3k(r)\ 

\j€Mn 


Then  measurement  z[l)  (  i  =  1,2  ,■■■  ,m(k))  is  validated  if  and  only  if 

rf,-«'(r)],[4'(r)]-1l4‘)-^W]<7  (17) 

where  7  is  an  appropriate  threshold.  The  volume  of  the  validation  region  with  the  threshold  7  is 

Vs(r)  :=  c„,7”'/2|Sj'(r)|l72  (18) 

where  nz  is  the  dimension  of  the  measurement  and  is  the  volume  of  the  unit  hypersphere  of 
this  dimension  (ci  —  2,C2  =  7r,C3  =  4tt/3,  etc.).  Choice  of  7  is  discussed  in  more  detail  in  [5,  Sec. 
2.3.2].  After  performing  the  validation  for  each  target  separately,  we  take  the  validation  region  for 
the  whole  target  set  as  the  union  of  separate  validation  regions  of  all  targets  .  As  mentioned  in 
[2]  and  [5],  "...  this  approach  is  adopted  in  order  to  have  the  PDF  of  each  false  measurement  the 
same,  i.e.,  uniformly  distributed  in  the  entire  validation  region”.  The  volume  of  validation  region 
for  the  whole  target  set  is  approximated  by 

Vk  =  Y. "Mr).  (19) 

r=l 

Step  1.4.  State  estimation  with  validated  measurements  (V)  €  Mn,  Vr  €  Tn)-  From 
among  all  the  raw  measurements  at  time  k,  i.e.,  Zk  :=  {2^,  z^k\  •  •  • ,  z^^},  define  the  set  of 
validated  measurement  for  sensor  1  at  time  k  as 

n==ai),yi2,.-.»r,)}  (20) 

where  m(k)  is  the  total  number  of  validated  measurement  at  time  k  and 

v[k  ■=  zki]  (21) 
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where  1  <  h  <  h  <  •  •  •  <  ljn(k)  <  m(k)  when  rfi(k)  0.  Note  that  all  targets  share  a  common 
validated  measurement  set  Y^.  , 

We  now  consider  joint  probabilistic  data  association  across  targets  following  [2]  and  [5].  Define 
the  validation  matrix 


ft  =  [o>jr]  i  —  !,•••  ,m{k),  r  =  0, 


(22) 


where  =1  if  the  measurement  i  lies  in  the  validation  gate  of  target  r,  else  it  is  zero.  A  joint 
association  event  ©a  is  represented  by  the  event  matrix 


ft(©fc)  =  P>ir(©*)]  *  =  1,*  ••»”»(*),  r  =  0,  •  •  • ,  JV 


where 


Uir(Qk) 


f  1  1 

1°  ‘ 


if  6ir{k)  C  ©A; 

otherwise. 


(23) 


(24) 


A  feasible  association  event  is  one  where  a  measurement  can  have  only  one  source 
f>ir(0fc)  =  l  Vi, 

r= 0 

and  where  at  most  one  measurement  can  originate  from  a  target 

m(k) 

^r(0fe)  :=  X)  wir(©fc)  <  1  for  r  =  l,---  ,N. 
i- 0 


(25) 


(26) 


The  above  joint  events  0£  are  mutually  exclusive  and  exhaustive. 

Following  the  definitions  in  [2]  and  [5],  define  the  binary  measurement  association  indicator 


N 


n(Qk)  :=  X^wir(0fc),  i  =  !,•••, m(k), 


(27) 


r= 1 


to  indicate  whether  the  validated  measurement  is  associated  with  a  target  in  event  0*.  Fur¬ 
thermore,  the  number  of  false  (unassociated)  measurements  in  event  ©a  is 


m(k) 

=  X  [1  -  T<(©fc)]. 

t=l 


(28) 


We  will  limit  our  discussion  to  nonparametric  JPDA  [2]  ,[5].  One  can  evaluate  the  likelihood  that 
the  target  r  is  in  model  j r  as 


A£r(r)  :=  p[Yk\Mlr (r), Zk~x]  =  X PiYk I ©fc , (r) , Zk~x]P{Qk | M3kr (r) , Zk~ 1 } 
=  X)p[nie*,^r(r),2:fc-1]p{©fc} 

Qk 


(29) 
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where,  as  in  [2, (9-31)],  the  irrelevant  conditioning  terms  have  been  omitted  in  the  last  line  of  (29) 
and  the  conditioning  on  fn(k)  is  implicit  in  the  event  0*,.  The  first  term  in  the  last  line  of  (29)  can 
be  written  as 

p[Yk\ek,Mt(r),Zk-']=±::  ±  £  ...  ± 

Jl=l  3r- 1  =  1  jr +1=1 

p[yfc|0fc,  Mi1  (1), . . . ,  Mi-1  (r  -  1),  Mi (r),  M*+1( r  4- 1),  •  •  • ,  M3kN (AT),  Z^1] 


xP{Mi(l),...,Mi-l(r-l),Mi+1(r  +  l),...,Mi”(N)\Gk,Mi(r),Zk-1}.  (30) 


The  second  term  (apriori  joint  association  probabilities)  in  the  last  line  of  (29)  turns  out  to  be  ([5, 
Sec.  6.2],  2,  Sec.  9.3]) 

p{efc}  =  ^^n(p^5(et)(1-^)1"Mefc)  (31) 

where  Pd  is  the  detection  probability  (assumed  to  be  the  same  for  all  targets)  and  e  >  0  is  a 
“diffuse”  prior  (for  nonpar ametric  modeling  of  clutter)  whose  exact  value  is  irrelevant.  We  assume 
that  the  states  of  the  targets  (including  the  modes)  conditioned  on  the  past  observations  are 
mutually  independent.  Then  the  first  term  on  the  right-side  of  (30)  can  be  written  as 

p[Yk\Gk,  Mi (1),  •  •  • ,  Mi~l  (r  -  1),  Mi (r),  M^ (r  +  1),  •  •  • ,  MjkN  (. N),Zk -1] 


m(k) 

II  P[y<£)\0iri(k),Mii(ri),Zk-1],  9iri(k)  C  0fc, 

i=l 


(32) 


where  the  conditional  PDF  (probability  density  function)  of  the  validated  measurement  y given 
its  origin  and  target  mode,  is  given  by 

f  if Ti(@k)  =  1, 

I  1/Vi  if  Ti(@k)  =  0 


P^k  \0*rt (k),  MJkT'  in),  zk  x] 


rfc-ll 


where 


M(x\ y,P )  :=  \2nP\  1/2exp[-^(a:  -  y)’P  x(a:  -  y)\. 


(33) 


(34) 


The  second  term  on  the  right-side  of  (30)  is  given  by 


P{M>' (1),  •  •  • ,  Ml’-' (r  -  1),  Mi’*' (r  + 1),  •  ■  • , M™ (N)\Sk, M>’ (r),  Z*-1} 

=  ft  r{M^W|et>M*(r),2t-1}=  n  pi  K(s)\Z“-l}=  0  /£»•  (35) 

s=l,s^r  s=l,s^r  s=l,s^r 
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The  probability  of  the  joint  association  event  0*  given  that  model  j  is  effective  for  target  r  from 
time  k  —  1  through  k  is 


P{ek\Mi(r),Zk~\Yk}  =  IplYkie^MlirlZ^PiQ^Mlir^Z^y 

=  \p[Yk\Qk,Mi{r),Zk-'}P{®k}  =:  $(r,0fc)  (36) 


where  the  first  term  can  be  calculated  from  (30)  and  (32)  -  (35),  the  second  term  from  (31),  and  c 
is  a  normalization  constant  such  that  Ylek  fik(ri  ©*)  =  1- 

The  following  updates  are  done  for  each  target  r  (r  €  Tn)-  Calculate  A kr(r)  (needed  in  Step 
1.5  later)  via  (29)-(35).  Define  the  target  and  mode-conditioned  innovations 


4(r>6k) 


yk }  -z{{r)  for  *  ==  1,  —  ,  m(k)  if  OiT{k)  C  0fe 
0  otherwise. 


(37) 


Using  ^.|fc_1(r)  (from  (12))  and  its  covariance  F^A._1(r)  (from  (13)),  one  computes  the  state  update 
k(r)  and  its  covariance  Pk\k{r )  according  to  the  standard  PDAF  [2], [5]: 

Kalman  gain:  W3k{r)  =  Pk\k^{r)H3k  (r)[S3k(r)]~l.  (38) 


State  estimate  update:  ar£|fe(r)  :=  E{xk(r)\M3k(r), Zk  \ T*}  =  ^ xj|fc(r,  Qk)/33k(r, @k) 


e* 

%(r,©fc)  :=  E{xk(r)\Ml(r),Zk-\Yk,ek}  =  ^(r)  +  t^'(r)^(r,©*). 
Covariance  of  54(0  =  H\k(r)  = 

ek  i=i 


(39) 

(40) 

(41) 


where 

Ai(Qk)  =  E  {xk{r)xk(r)\MJk(r),@k,Zk-\Yk}  =  x3klk{r,ek)x{]k(r,@k)  +  P*\k(r,Qk),  (42) 


•*4(r’®fc)  — 


H\k-i  (r)-Wl(r)S{(r)Wi'(r) 

Pk\k-i(r) 


if  9ir(k)  C  ©fc  for  some  1  <  t  <  m(k) 
otherwise, 


(43) 


A2(©fc)  =  “^(r.0t)®(n(r),  A3(0*)  =  -%Mx{|k(r>  0*)>  ^(©fc)  =  K\k{f)^k\k{r)-  (44) 


Step  1.5.  Update  of  mode  probabilities  (Vj  G  Mn,  Vr  G  Tn): 

d(r)  ==  P[M{{r)\Zk)  =  P[M3k{r)\Zk-l)p[Yk\M3k{r\  Zk~l\  =  -/k~{r)K 3k(r)  (45) 

where  c  is  a  normalization  constant  such  that  i  /ii(r)  =  !• 
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Step  1.6.  Combination  of  the  mode-conditioned  estimates  (Vr  €  Tn)‘  The  final  state 
estimate  update  at  time  k  is  given  by 

**!*(»•)  =  ]C*fc|fc(rK(r)  (46) 

i= 1 

and  its  covariance  is  given  by 

pk\k(r)  =  it  { P3k\k{r )  +  [xi\k{r)  -  *fc|fc(r)][2j^(r)  -  £*|*(r)]'}  /*J(r).  (47) 

3= 1 

IV.2  Second  Scan  Steps 

Here  we  update  to  scan  k  +  1,  given  data  up  to  time  k  +  1,  with  a  sliding  scan  window  of  size  two, 
scans  {k,  k  +  1}. 

Step  2.1.  Interaction  —  mixing  of  the  estimate  from  the  previous  time  (Vj  G  Fin, 
Vr  G  Tat): 

Ml(r,&k)  :=  P{Mt(r)\Z\ek}  =  ^p[Yk\Mi{r),Zk-\@k)P{Mi{r)\Zk-\ek} 

=  cPl(r,©k)n{~(r).  (48) 

n 

predicted  mode  probability:  /ii+i (r,  Ok)  :=  P  {M{+l{r)\Zk ,  ©*}  =  ^py-/4(r,  ©*)•  (49) 

i= 1 

mixing  probability:  ^j(r,  ©*)  :=  P{Ml(r)\M3k+1(r),  Zk,  0fc}  =  Pijpi(r,  0*)//4+i(r>  ©fc).(50) 
mixed  estimate:  ^(r,©*)  :=  E{xk{r)\Ml+1(r),Zk,Qk}  =  J2Bk\k(r>®kWlj(r,®k)-  (51) 

i- 1 

covariance  of  the  mixed  estimate: 

P%k(r,Vt)  EfeM  -  S^(r,e4)]j*»(r)  -  *J{t(r,et)]'|i4+,(r),*,eU 
=  4  eetn(r,e*)  -  £&(>-,  eoMi/blr,©*)  -s&fr.eoiV'tr.e*).  (52) 

i—l 

Step  2.2.  Predicted  state  (Vj  G  Mn,  Vr  G  Tn)\ 

State  prediction:  xj[+i|fc(r’ 0*)  :=  Eixk+i(r)\M3k+1(r), Zk, 0fe}  =  F3kx^k{r, ©*)•  (53) 

State  prediction  error  covariance: 

Pfc+i|fc(r>0*)  =  p([*M-i(r)  -  ©*)][®*+i(»*) -&k+\\k(r'ek)]'\M3k+l{r),Zk,Qk} 

=  HP°k!k(^k)Fl'+GiQiG{.  (54) 
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Using  (2)  and  (53),  the  mode-conditioned  predicted  measurement  of  target  r  is 

Sk,(r,e»):=V(3k1|t(  r,e4».  (55) 

Using  the  linearized  version  (3),  the  covariance  of  the  mode-conditioned  residual  4+ \(r,®k)  '■= 

4+i  -  **+1  (r>  ©fc) is  eiven  by 

Sl+1(r,ek)  :=  ^{4+i(r>efe)4+i(r>0fe)'lAffe+i(r)>'2fc>0fc} 

=  <1(r,0fe)P^1|fc(r,©fc)<1(r,  ©*)  +  <!  (56) 

where  Hl+l(r,Qk)  is  the  first  order  derivative  (Jacobian  matrix)  of  hJ(.)  at  4+i|fc(r’  ©*)• 

Step  2.3.  Measurement  validation:  For  target  r,  the  validation  region  is  taken  to  be  the 
same  for  all  models  and  0fe’s,  i.e.,  as  the  largest  of  them.  Let 

(jr,©fc)  “  arg  (^max^^  |S£+1(r,0fc)|)  •  (57) 

Then  measurement  4+i  (  *  =  1, 2,  •  •  •  ,m{k  -t- 1))  is  validated  if  and  only  if 

[4+1  -  4+1  (r>  ©*)]  W+i  (r>  ©fc)]_1  [4+1  -  4r+i(r’  ©a)]  <  7  (58) 

where  7  is  an  appropriate  threshold.  The  volume  of  the  validation  region  with  the  threshold  7  is 

Vk+i(r)  :=  (^zjnz/2\Slr+1(r,ek)\1/2  (59) 

where  nz  is  the  dimension  of  the  measurement  and  criz  is  the  volume  of  the  unit  hypersphere  of 
this  dimension.  The  volume  of  validation  region  for  the  whole  target  set  is  approximated  by 

V*+I  =  EWr).  (60) 


Step  2.4.  State  estimation  with  validated  measurements  (Vj  €  _Mn,  Vr  6  Tn):  From 
among  all  the  raw  measurements  at  time  k- f  1,  i.e.,  Zk+i  :=  {4+i>  4+l> ' "  >  4+ifc+1^}>  define  the 
set  of  validated  measurement  for  sensor  1  at  time  k  +  1  as 


(ro(fe+l))i 


(61) 


where  m(k- f  1)  is  the  total  number  of  validated  measurement  at  time  fc  +  1.  and 

4+1  :=  4+1  (62) 

where  1  <  /i  <  h  <  ••  •  <  lfn(k+i)  <  wi(A  +  1)  when  m(k  +  .1)  ^  0.  Note  that  all  targets  share  a 
common  validated  measurement  set  Yk+i- 
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We  now  consider  joint  probabilistic  data  association  across  targets  as  in  Sec.  IV.  1.  Define  the 
validation  matrix 


fi  =  [wir]  *  =  1,  •  ••  ,m(k  +  1),  r  =  0,---,N 


(63) 


where  uj{r  =1  if  the  measurement  i  lies  in  the  validation  gate  of  target  r,  else  it  is  zero.  A  joint 
association  event  0*+i  is  represented  by  the  event  matrix 

^(0/fc+i)  =  [tiiT{®k+i)]  *  =  l,---,m(fc  +  1),  r  =  0,  •  •  • ,  TV  (64) 


where 


Wir(©fc+l)  = 


1  if  0ir(&  +  l)  C  0fc+i 
0  otherwise. 


(65) 


A  feasible  association  event  is  one  where  a  measurement  can  have  only  one  source 
N 

]^£tr(©fc+i)  =  1  Vi,  (66) 

r=0 

and  where  at  most  one  measurement  can  originate  from  a  target 
m(fc+ 1) 

<*r(0fc+i)  :=  S*r(©fc+i)  <  1  for  r  =  1,  ■  •  • ,  JV.  (67) 

»=o 

The  above  joint  events  ©*+ 1  are  mutually  exclusive  and  exhaustive. 

As  in  Sec.  IV.  1,  define  the  binary  measurement  association  indicator 
N 

Ti(©jfc+i)  :=  ^£Dir(0fe+i),  *  =  !,•■•  ,m(k  + 1),  (68) 

r= 1 

to  indicate  whether  the  validated  measurement  is  associated  with  a  target  in  event  0*.+1. 
Furthermore,  the  number  of  false  (unassociated)  measurements  in  event  ©*4. i  is 
m(fc+l) 

<f>(ek+i)=  [1  -  Tf(©*+1}].  (69) 

t=i 

We  will  limit  our  discussion  to  nonparametric  JPDA  [2], [5].  One  can  evaluate  the  likelihood  that 
the  target  r  is  in  model  jr  as 


K+dr)  :=pK+i|M>  ,(r),Z‘] 


=  E  E  p[n+i|0fc,©fe+i,Mf+1(r),^]P{0fc+1}P{0fc|M>  x(r),2fc}.  (70) 

©*  ©/fc+i 
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The  first  term  in  the  last  line  of  (70)  can  be  written  as 

p[yfc+1|e,,efc+1,M^1(r),zfc]=  E  E  •••  E 

jl=l  jr-l=lir+l=l  iiV=l 

pWh-iI©*,  ©*+i,-Mfc+i(l),-  •  • , Mjfc1  (r  -  1), Mf+1(r),  Mjfc1 (r  +  1),  •  •  • ,  M^N),  Zk] 

xP{M*_x (1),  •  •  • ,  Affc+i1  (r  -  1), M/;Y  (r  +  1),  •  •  • ,  ^(JV)!©*,  0*+1,  M>+1  (r),  2*}.  (71) 

The  second  term  (apriori  joint  association  probabilities)  in  the  last  line  of  (70)  turns  out  to  be  ([5, 
Sec.  6.2],  2,  Sec.  9.3]) 

p{©*+1}  =  i^iyf  mpo)s^Hi  -  (72) 

where  Pd  is  the  detection  probability  (assumed  to  be  the  same  for  till  targets)  and  e  >  0  is  a 
“diffuse”  prior  (for  nonparametric  modeling  of  clutter)  whose  exact  value  is  irrelevant.  The  third 
term  in  the  last  line  of  (70)  is  given  by 

P{Qk\Mt+1(r),Zk}  =  c'P{M£+l(r)\ek,Zk}P{&k\Zk}  =  ctf-.ir, Qk)P{ek}p\Yk\ek, -Z*'1]  (73) 
where 


p[Yk\Qk,Zk-1}  =  J2p[Yk\ek,Mi(r),Zk-1]4-(r).  (74) 

3 

We  assume  that  the  states  of  the  targets  (including  the  modes)  conditioned  on  the  past  observations 
are  mutually  independent.  Then  the  first  term  on  the  right-side  of  (71)  can  be  written  as 


P[Yk+i I©*,  ©fc+i,  m£+1(  1), . . . ,  M{^  (r  -  1),  Af* x(r),  « (r  +  1),  •  •  • ,  M^N),  Zk] 

m(/c+ 1) 

*  II  PlPfc+il^r4(fc  +  l),©fc,M^1(r<),^],  ftrj(Hl)c0Wl 

t=l 


(75) 


where  the  conditional  PDF  of  the  validated  measurement  given  its  origin  and  target  mode,  is 
given  by 


>< 


rk]  _ 


M(vk+v  K+1  (ri>  ©fe).  ^+i(n,©fc))  if  Ti(0fc+i)  =  l, 

if  Ti(0ifc+i)  =  0. 


l/Vfc+1 

The  second  term  on  the  right-side  of  (71)  is  given  by 

=  n  ,(«)!©!+,, et.M*  ,^,2*}=  n 

5=1,  s^r  5=1,  s^r 


(76) 
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N 


=  n  /4s+i(s’0*)- 

s=l,57^r 


(77) 


The  probability  of  the  joint  association  events  0fc+i  and  @k  given  that  model  j  is  effective  for 
target  r  from  time  k  through  k  +  1  is 

P{Gk+u  ek\Mi+1(r),  Zk,  Yk+1}  =  ip[n+1|0A+1,  ©*,  Ml+1(r),  Zk]P{ek+1}P{ek\Mi+1(r),  Zk) 


-■P{+ i(r>  ©*+!>©*) 


(78) 


where  the  first  term  can  be  calculated  from  (71)  and  (75)  -  (77),  the  second  term  from  (72),  the  third 
term  from  (73),  and  c  is  a  normalization  constant  such  that  J2ek+1  X)©fc  i®jt+i(r>  ©fc+ii  ©fc)  =  1- 
The  following  updates  Eire  done  for  each  target  r  (r  €  Tn)-  Calculate  A£r+1(r)  (needed  in  Step 
2.5  later)  via  (70)-(77).  Define  the  target  and  mode-conditioned  innovations 


4+i(‘ 


',©*+!,©*):=<  V‘ 

\  0 


_  /  yl+i  “  ^+i(r>  ©fc)  for  *  =  m{k  +  1)  if  6ir(k  + 1)  C  0fc+i 


otherwise. 


(79) 


Using  2*+1|fc(r,  @k)  (from  (53))  and  its  covariance  Pk+uk(r,  Qk)  (from  (54)),  one  computes  the  state 
update  xj[+1|fe+1(r)  and  its  covariance  P^+1|fc+1(r)  as  follows. 


Kalman  gain:  Wj^r,©*)  =  P£+yk{r,Gk)H{+1(r,ek)[Si+1(r,  @k)]  \ 
State  estimate  update:  ar£+1|fc+1  (r)  :=  E{xk+ 1 (r) | M3k+1  (r),  Zk,  Yk+i} 

~  El  EZ  a4+iifc+i(7’’  ©fc+ii  ©fc)/^fc+i(r)  ©fc+i)  ©A:) 

©fc+i  @jb 


(80) 


(81) 


4+i|*+i(r>0fc+i,©fc)  =^+i|fc(^,©fc)  +  Wfc+i(r,0*)i4+i(r-,©fe+i,0jb).  (82) 

4 

Covariance  of  4+i|Jfc+i(r)  :  ^fc+ilJfc+i (r)  =  E  EE M@k+ 1, ©*)^+1(r,  ©*+x, 0fc)  (83) 

©fc+l  ©Jk  <=1 

where 

■Ai(©fc+i,©*)  =  E{xk+i{r)x'k+1(r)\M3k+l{r),@k+1,Gk,Zk,Yk+i} 

= ®i+iifc+i  (r,  ©*+i,0fc)  +  P^+1|fc+1(r,  ©fc+i,  ©*),  (84) 

■^fc+i|fc+i(r’  ©fc+i»©fe) 

_  I  i^i|*(r,0*)  -  W3k+1(r,@k)Sl+1(r,ek)Wl'+1(r,ek)  if  0ir{k  + 1)  C  0fc+i,  1  <i<  m{k  +  1) 


JW'.0*) 


otherwise, 


(85) 
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-^(©fc+li  ®k)  —  a4+l|fc+l(r’ ®fc)*Efe+l|jt+l(r)> 

(86) 

A3(0fe+i,  @fc)  =  -^+i|fc+i(r)^,+1|fc+1(r,0*+i,0fc), 

(87) 

Ai(0fc+i,©fc)  =  4+i|jfe+i(r)4+i|fe+i(r)- 

(88) 

Step  2.5.  Update  of  mode  probabilities  (Vj  E  Mn,  Vr  E  Tn)'- 

«Ui(r)  :=  i’W+,(r)|2‘+1]  =  PK+1M|2*Wn+l|M'+l(r),2‘]  =  ^£i(r)4+,(r)  (89) 

where  c  is  a  normalization  constant  such  that  /4+1  (f)  —  1  and 

viii(r)  =  (90) 

e* 

Step  2.6.  Combination  of  the  mode-conditioned  estimates  (Vr  €  Tn)'  The  final  state 

estimate  update  at  time  k  +  1  is  given  by 
n 

®fc+i|*+i(r)  =  ESfc+i|*+i(rH+i(r)  (91) 

j=l 

and  its  covariance  is  given  by 

Pk+i\k+i(r) 

n 

=  E  {pk+i\k+i(r)  +  [**+i|*+i(r)  -  ®fe+i|*+i (*-)][®fc+i|fc+i (r)  -  **+i|ib+i(r)]'}  /4+i(r).  (92) 

j= i 

V  Simulation  Example 

We  now  consider  tracking  three  maneuvering  targets  in  clutter.  We  carry  out  state  estimation  for 
each  target  using  IMM  multiscan  JPDA  with  a  scan  window  size  of  two  and  compare  our  results 
with  single  scan  IMM/ JPDA  algorithm  of  [8]. 

The  True  Trajectories:  Target  1  starts  at  location  [10500  1740  40]  in  Cartesian  coordinates 
in  meters.  The  initial  velocity  is  [-140  299.9  0]  in  m/s.  Target  stays  at  constant  altitude  with  a 
constant  speed  of  331  m/s.  Its  trajectory  is  a  straight  line  with  constant  velocity  between  0  and 
15  sec.,  a  coordinated  turn  of  -0.32  rad/s  with  a  constant  acceleration  of  109  m/s2  between  15  and 
25  s,  and  a  straight  line  with  a  constant  velocity  between  25  and  35  s.  Target  2  starts  at  location 
[9800  1960  40]  in  Cartesian  coordinates  in  meters.  The  initial  velocity  is  [0  299  0]  in  m/s.  The 
target  stays  at  a  constant  altitude  with  a  constant  speed  of  299  m/s.  Its  trajectory  is  a  straight 
line  with  constant  velocity  between  0  and  15  sec.,  a  coordinated  turn  of  0.32  rad/s  with  a  constant 
acceleration  of  94  m/s2  between  15  and  25  s,  and  a  straight  line  with  a  constant  velocity  between 
25  and  35  s.  Target  3  starts  at  location  [9200  1740  40]  in  Cartesian  coordinates  in  meters.  The 
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initial  velocity  is  [0  299  0]  in  m/s.  The  target  stays  at  a  constant  altitude  with  a  constant  speed  of 
299  m/s.  It’s  trajectory  is  a  straight  line  with  a  constant  velocity  between  0  and  35  sec. 

The  Target  Motion  Models:  The  motion  models  are  identical  for  all  three  targets.  In  each 
mode  target  dynamics  are  modeled  in  cartesian  coordinates  as  xb  =  Fxk- 1  +  Gvk- \  where  state 
of  the  target  is  position,  velocity  and  acceleration  in  each  of  the  three  Cartesian  coordinates  ( x ,  y 
and  z).  Thus  Xk  is  of  dimension  9  (nx=9).  Three  models  are  considered  in  the  following  discussion. 
The  system  matrices  F  and  G  are  defined  as 
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•  Model  1:  nearly  constant  velocity  model  with  zero  mean  perturbation  in  acceleration. 
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where  T  is  the  sampling  period.  The  standard  deviation  of  the  process  noise  of  M1  is  5m/s2. 
•  Model  2:  Wiener  process  acceleration  (nearly  constant  acceleration  motion) 
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The  standard  deviation  of  the  process  noise  of  M2  is  7.5m/s2. 

•  Model  3:  Wiener  process  acceleration  (model  with  large  acceleration  increments,  for  the  onset 
and  termination  of  maneuvers).  Here  Fjf  =  Ffc2  and  G\-Gl.  The  standard  deviation  of  the 
process  noise  of  M 2  is  40m/s2. 

The  initial  model  probabilities  for  three  targets  are  identical:  =  0.8,  yl  =  0.1  and  /i$  =  0.1.  The 

mode  switching  probability  matrix  for  three  targets  is  also  identical  and  is  given  by 


Pll 

P12 

P13 

'  0.8 

0.1 

O 

P21 

P22 

P23 

0.1 

0.8 

0.1 

_  P31 

PZ2 

P33 

0.1 

0.1 

0.8 

The  Sensor:  A  single  sensor  (radar)  is  used  to  obtain  the  measurements.  The  measurements 
are  range,  and  azimuth  and  elevation  angles  for  the  radar).  The  range,  azimuth  and  elevation 
transformations,  respectively,  are  given  by 

r  =  (x2  +  y2  +  22)1/2,  a  —  tan_1(y/z),  e  —  tan_1[z/(x2  +  y2)1/2]. 
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The  measurement  noise  uij,  is  assumed  to  be  zero-mean  white  Gaussian  with  known  covariance 
matrix  R  =  diag[400m2,  49mrad2,  4mrad2].  The  sensor  is  assumed  to  be  located  at  the  origin  of 
the  coordinate  system.  The  sampling  interval  was  T  =  Is  and  it  was  assumed  that  the  probability 
of  detection  Pd  =  0.997. 

The  Clutter:  For  generating  false  measurements  in  simulations,  the  clutter  was  assumed  to 
be  Poisson  distributed  with  expected  number  of  A  =  0.1  /m  rad2.  These  statistics  were  used  for 
generating  the  clutter  in  all  simulations.  However,  a  nonparametric  clutter  model  was  used  for 
implementing  all  the  algorithms  for  target  tracking. 

Other  Parameters:  The  gates  for  setting  up  the  validation  regions  for  the  sensor  were  based  on 
the  threshold  7  =  16.  With  the  measurement  vector  of  dimension  3,  this  leads  to  a  gate  probability 
PG  =  0.9989  (see  p.  96  of  [5]). 

Simulation  Results:  The  results  were  obtained  from  30  Monte  Carlo  runs.  Fig.  1  shows  the 
true  trajectories  of  the  three  targets  and  the  distances  among  targets  as  a  function  of  time.  Fig.  2 
shows  the  RMSE  (root  mean-square  error)  for  the  filtered  position  estimates  for  the  three  targets 
as  a  function  of  time.  It  is  seen  from  Fig.  2  that  the  multiscan  approach  does  provide  a  significant 
improvement  over  the  single  scan  approach. 

VI  Conclusions 

We  investigated  the  problem  of  tracking  multiple  maneuvering  targets  in  the  presence  of  clutter 
using  switching  multiple  target  motion  models.  A  novel  suboptimal  filtering  algorithm  was  de¬ 
veloped  by  applying  the  basic  interacting  multiple  model  (IMM)  approach  and  multiscan  joint 
probability  data  association  (JPDA)  technique.  Past  work  (see  [14])  on  this  problem  is  restricted 
to  non-maneuvering  targets. 

The  algorithm  was  illustrated  via  a  simulation  example  involving  tracking  of  three  maneuvering 
targets  and  a  multiscan  data  window  of  length  two.  The  simulation  example  shows  a  significant 
improvement  in  target  position  estimate  by  the  proposed  IMM  multiscan  JPDA  (with  a  scan  window 
size  of  two)  compared  to  the  results  of  the  single  scan  IMM/ JPDA  algorithm  of  [8]. 
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Figure  1:  Trajectories  of  the  three  maneuvering  targets  (read  left  to  right,  top  to  bottom),  (a) 
Position  in  xy  plane,  (b)  distance  between  target  pairs,  (c)  x  and  ^velocities,  (d)  acceleration 
magnitudes. 
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Figure  2:  Root  mean-square  error  (RMSE)  in  position  using  single  scan  IMM/JPDAF  [8]  and  the 
proposed  multiscan  (window  size  2  scans)  IMM/Mscan-JPDAF  algorithms. 
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